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Abstract 


A discussion  of  absorption  phenomena  in  the  infrared  (2.5  - 
15pm)  spectral  region  is  presented.  Mechanisms  surveyed  in- 
clude the  attenuation  of  electromagnetic  waves  by  dielectrics 
and  metals;  molecular  absorption  processes;  lattice  vibrations 
in  bulk  and  powdered  dielectrics;  band  gap,  free  carrier  and 
multiple  phonon  absorption  in  semiconductors;  absorption  in  bulk 
and  powdered  metals;  and  cibsorption  by  polymeric  one-dimensional 
conductors.  Also  presented  is  a review  of  quantitative  methods 
of  determining  the  optical  constants  of  solids,  liquids  and 
gases  cuid  their  possible  sources  of  error.  The  results  of  quan- 
titative cibsorption  measurements  on  tri-n-butyl  phosphate, 
polymethylmethacrylate,  cellulose,  polyvinyl  alcohol  and  urea- 
formaldehyde  are  included. 
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SURVEY  OF  ABSORPTION  PHENOMENA  IN  THE  NEAR  INFRARED  REGION 


I . 


A.  Attenuation  of  Electromagnetic  Waves  by  a Dielectric 
Medium:  A Simplified  Dispersion  Theory 

We  would  like  to  describe,  by  the  use  of  a simple  model 
the  absorption  of  electromagnetic  radiation  by  non-conducting 
matter.  We  consider  the  material  as  consisting  of  a collection 
of  N identical  one-dimensional  oscillators  of  mass  m,  charge  e 
and  elastic  constant  q,  per  unit  volume.  Free  oscillations  are 
described  by 

mii  + mal^*u  = 0 (1) 

in  which  u is  the  displacement  of  the  particle  from  its  equi- 
librium position,  m is  the  mass  and  = q/m  is  the  reso- 
nant frequency.  (The  charge  does  not  appear  at  this  point.) 

If  the  oscillator  is  coupled  to  a loss  mechanism,  a damping 
term  must  be  added  to  the  above  equation  to  account  for  the 
energy  loss.  The  resulting  damped  oscillator  equation  is 

mu  + myu  + mw^u  = 0 (2) 

in  which  y is  the  damping  constant. 

The  resonauice  frequency  is  different  from  u)^.  The  inter- 
action of  the  incident  electromagnetic  field  with  the  oscillator 
is  assumed  to  be  as  follows:  The  periodic  electric  field  of  the 
incident  wave  exerts  a force  on  the  damped  oscillator  and  drives 
it  at  a frequency  u).  We  describe  this  motion  with  the  forced 
oscillator  equation 

mii  + myu  + mw^^u  =*  Force  = (3) 

where  is  the  maximum  amplitude  of  the  incident  electromagnetic 
wave  oscillating  in  the  u direction  at  the  site  of  the  oscilla- 
tor . 
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We  can  solve  the  above  equation  by  substituting  a trial 
solution  of  the  form 


u(t) 


u e 
0 


-iu)t 


(4) 


into  equation  (3) . The  result  is 

-iwt 


u(t)  = 


eE  /m 
o 


e/m 


- iyiD 


oj  ^ - 0)*  - ivo) 
0 


E(t)  (5) 


Thus,  the  amplitude  u depends  on  the  parameters  of  the  oscil- 
lator (e,m,a)^,Y»)  as  well  as  on  the  strength  of  the  incident 
electromagnetic  field. 

So  far  we  have  considered  just  one  microscopic  process. 

The  description  of  electromagnetic  effects  by  Maxwell's  equations 
uses  macroscopic  quantities  and  material  constants  such  as  the 
dielectric  constant  e and  permeability  y.  The  connection  be- 
tween the  microscopic  parameters  of  the  oscillators  and  e and  y 
is  of  interest  here. 

Since  the  amplitude  u,  derived  above,  is  the  displacement 
of  one  oscillator  of  mass  m and  charge  e,  eu  is  the  induced 
dipole  moment.  The  contribution  of  all  N induced  dipole  moments 
per  unit  volume  in  the  material  is  the  polarization 

P = Neu  (6) 

Introducing  u from  equation  (5)  above  gives 


m u)  ^ - 0)*  - iojY 
o 


It  should  be  mentioned  that  in  general  E is  a function  of  both, 
position  (x)  and  time  (t) . 

We  now  wish  to  solve  the  wave  equation  in  the  medium  mod- 
elled by  the  damped  oscillators.  From  equation  (A- 10)  in  appen- 
dix A,  for  a one-dimensional  dielectric  medium  with  no  conduc- 
tivity, the  appropriate  wave  equation  is 


3 


3 3 . 3 

= c u + u 

3 ° ° 3 ° 3 


(8) 


where  we  have  used  the  identity 


The  second  term  on  the  right  hand  side  of  equation  (8)  is  called 
a source  term.  It  accounts  for  the  effect  of  polarization 
charges  in  the  medium  as  discussed  in  appendix  A.  We  now  sub- 
stitute for  P and  E in  equation  (8)  displaying  explicitly  the 
spatial  and  temporal  dependence  of  the  electric  field,  i.e. 


where  K 
suit  is 


or 


E(x,t)  = E^e 


i {Kx-(jL)t) 


(10) 


27r 


is  the  magnitude  of  the  wave  vector.  The  re- 


-K*  = - - y 

0 0 o 


Ne‘ 


(u)^  -(*)*-  iojy 


(i) 


Ne  */me 


1 + 


(1) 


* - 0)^  - iooYj 


(11) 


again  using  equation  (9).  We  see  that  is  complex.  Since  K is 
related  to  the  index  of  refraction  by  K = n,  we  introduce 

C 

the  complex  index  of  refraction  n*=n+i)c. 

Then 

K*^  ^ (n  + ik)  ^ ^ + 2in)c 

U)2 


From  equations  (11)  and  (12)  we  obtain  the  real  and  imaginary 
parts  of  K**  as 

Ne*  (to  * - to*) 

n*  - )c*  = 1 + 2 (13) 

e^m  [(Wg^  - to*)*  + to*Y^] 
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and 


2nk  = 


Ne‘ 


e m 
o 


u)Y 


- 0)^2  + 


(14) 


We  call  n the  refraction  coefficient  and  k the  extinction  co- 
efficient. 

We  have  assumed  above  that  the  driving  electrical  field  at 
the  site  of  the  oscillator  is  equal  to  the  field  of  the  electro- 
magnetic wave  outside  of  the  material.  This  assumption  implies 
that  the  density  of  oscillators  in  the  material  is  low,  somewhat 
as  it  would  be  in  a gas.  Conversely,  for  a dense  distribution 
of  oscillators  we  would  have  to  add  the  influence  of  the  sur- 
rounding material,  the  so-called  Lorentz  contribution.  This 
latter  case  will  not  be  considered  here. 

For  such  a low  density  medium  the  refractive  index  must  be 
close  to  1,  the  value  in  free  space.  As  a consequence  k is 
small  compared  to  1 and  we  have  for  n and  k,  using 


/I  + X = 1 + i 


and 


n = 1 + 


k = 


Ne‘ 


2e  m 
0 


Ne‘ 


2e  m 
o 


- 0)2) 


(O)  2 - U)2  ) ^ + U)2y  ‘ 


(u). 


2 _ ^2)2  ...  ...2.,2 


(15) 


(16) 


+ 0)2y' 


It  is  instructive  to  plot  n(o))  and  k (o))  and  to  consider  three 
specific  frequencies,  o)  = 0,  u)  = u)^  and  o)  = ■».  To  facilitate 
this  discussion,  we  introduce  the  frequency-like  term  0)p=Ne2/e^m. 
Then  equations  (15)  and  (16)  become  respectively 


n (o))  = 1 + Ty  0)^2 
2 P 


0)  2 - oj* 

0 


(0)^2  - (1)2  )2  + y*0)2 


(17) 
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Consider  first  the  refractive  index,  n(a))  plotted  in  Fig- 
ure 1 . 


0) 

0 


Fig.  1.  Frequency  dependence  of  refractive  index  for  a 
dielectric  described  by  a simplified  dispersion 
mode 1 . 

For  u)  = 0,  and  y = 0 (corresponding  to  a lossless  medium, 
n(0)  = 1 + *5  (ojp/o)^)  ^ = 1 + (Ne^/2e^q).  Numerically,  n(0) 

should  correspond  to  expected  static  values  for  low  density 
media,  i.e.,  n(0)~l.  To  see  this,  we  choose  reasonable  vai-fis 

of  N and  q and  evaluate  (Wp/u)^)*,  i.e.. 


For  (!)-*•«> 


Thus  we  see  that  n(0)  = 1 as  expected. 


, n(«>)+l 


which  is  the  expected  high  frequency  limit  (for  example,  in  the 
x-ray  region).  For  u)  = (and  y = 0)  , n((i))  is  singular.  This 
situation  is  never  realized  for  real  systems,  however,  since  some 
loss  mechanism  is  always  present.  The  effect  of  the  non-zero 


damping  term  on  n(a))  is  shown  as  a dashed  line  in  Figure  1.  For 
this  latter  case,  when  w ~ n(o3)  > 1 and  is  increasing  as  w 

increases.  This  is  normal  dispersion  behavior  and  is  illustrated 
when  white  light  passes  through  a glass  prism  with  violet  light 
bending  more  than  red  light.  Near  resonance  (w  = 

n (w)  decreases  as  oj  increases  giving  rise  to  anomalous  dispersion. 
This  effect  would  be  observed  as  the  opposite  as  that  just  des- 


cribed, i.e.  red  light  is  bent  more  than  blue  in  traversing  the 


glass  prism. 

We  now  consider  the  extinction  coefficient  k.((i)).  From 


equation  (18),  k(0)  = 0 and  k (“)  -*•  0.  At  resonance  k(aj^)  has  a 
maxima,  the  width  depending  upon  y»  the  damping  term.  Clearly 
for  y = 0,  k(w)  = 0 for  all  w.  k (to)  is  plotted  in  Figure  2. 


Fig.  2.  Frequency  dependence  of  extinction  coefficient  for 
a dielectric  described  by  a simplified  dispersion 
model . 
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It  should  be  emphasized  that  the  optical  behavior  just 
described  was  based  upon  a very  elementary  model  composed  of  N 
identical  oscillators.  Real  materials  require  more  realistic 
[TKjdels  to  accurately  predict  or  describe  observed  measurements. 
As  an  example  of  a first  attempt  to  extend  the  above  model,  we 
could  relax  the  restriction  of  identical  oscillators  and  assume 
that  there  are  a certain  fraction  f^  of  particles  having  reso- 

(ij  . = / q . /m. 

rk  T ^ 1 ' 


nant  frequencies 
becomes 


OD 


Then,  for  example,  equation  (11) 


K 


*2  = 


(O' 


Ne  ‘ 


e m 
o 


- 0) 


- iy-iii) 


-] 


(11') 


with  corresponding  expressions  for  n^-k^  and  2nk  following,  as 
before. 


Absorption  Coefficient  a 

We  now  consider  the  effect  of  a complex  wavenumber  K*  on 
the  propagation  of  a plane  electromagnetic  wave  in  a dielectric 
medium.  Equation  (10)  is,  using  K*=  ^ ^ 


i [p(n+ik) x-wt] 

E(x,t)  = E^e 


- ^x  i (^x-wt) 

= e ^ E^e  (19) 

Since  the  intensity  I of  the  radiation  is  the  square  of  the 
electric  field,  we  get 

-2(p  k)x 

I = )e|^  = EE*  = |e^|*  e ^ (20) 

Letting  = I | * and 

a * 2a  = 2 ^ k (21) 
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we  qet  Beer ' s Law 


I (x) 


I e 


o 


-ax 


(22) 


Thus  a tells  us  how  much  a layer  of  medium  x units  thick 
attenuates  a beam  of  radiation  of  intensity  incident  upon 

It.  In  particular,  a , the  absorption  coefficient,  measured 
in  units  of  reciprocal  centimeters,  is  such  that  in  a distance 
(a)"*  the  intensity  falls  from  to  I (=  I^/e) . Experimentally 
a is  measured  directly.  (See  discussion  below.) 

Using  equation  (21)  and  the  well-known  relations  connecting 
the  various  wave  parameters. 


0)  = 2'nf  , C = Af 


(23) 


(C  = speed  of  light  in  vacuo  = 3 x 10®  m/s ; X = wavelength  of 
radiation)  we  see  that 


and 


(24) 


a is  the  quantity  usually  called  the  absorption  coefficient. 
Equation  (24)  shows  that  a is  sufficient  to  determine  the  ex- 
tinction coefficient  k.  Note  that  in  entering  a medium  of  re- 
fractive index  n,  a plane  wave  (1)  slows  down  by  a factor  of 
1/n  (since  V = C/n) , (2)  the  wavelength  decreases  by  a factor  of 
1/n,  (^j^=  ^0  ~ space  wavelength)  , but  (3)  the  fre- 

quency remains  constant,  that  is 


U)  = 2Tlf 


(C/n) 


2ttC 

X 

0 


0) 


(25) 


B.  Attenuation  of  Electromagnetic  Waves  by  a Conducting 
Medium 


Let  us  compare  a conducting  medium  with  the  dielectric 
medium  previously  considered.  In  a dielectric  medium,  the 
bound  charges,  represented  as  oscillators,  form  induced  electric 
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¥ 


dipole  when  exposed  to  incident  electromagnetic  wave.  In  a 
conductor,  however,  this  polarization  mechanism  is  missing  since 
the  electrons  are  not  bound  by  a restoring  force,  but  are  con- 
sidered to  be  almost  free,  except  at  the  boundaries.  The  inci- 
dent electromagnetic  wave  now  does  not  produce  oscillations  of 
t.ound  electrons,  but  does  produce  microscopic  currents.  We 
therefore  expect,  in  the  wave  equation,  a current  density  re- 
lated term.  From  equation  (A-10)  of  appendix  A we  obtain,  for 
a one-dimensional  conducting  medium,  a wave  equation  of  the  form 


3 

a x^ 


y e 
o o 


3 

a t" 


+ 


li 

a t 


(26) 


where  j is  the  current  density  and  is  related  to  the  velocity 
of  the  electrons  by  j = Nev  or  to  the  applied  instantaneous 
electric  field  by  j = oE , o being  the  electrical  conductivity. 
Again  our  objective  is  to  correlate  the  conductors  microscopic 
parameters  with  the  measured  optical  "constants". 

If  we  try  to  apply  the  damped  oscillator  equation  (with  no 
restoring  force)  we  would  have 

rnii  + nryu  = eE^E  (27) 

The  general  solution  of  this  equation  is  the  sum  of  the  solution 
of  the  homogeneous  and  inhomogeneous  equations.  For  the  homo- 
geneous equation  we  have 

mvi  + myu  = 0 (28) 


We  try  as  a solution 


u = e 


t 

T 


and  obtain  Y = The  damping  constant,  y,  is  inversely  pro- 

portional to  a decay  time  t measured  from  an  arbitrary  zero 
point  of  time.  This  time  is  short  (typically  ~ 10  ' * sec)  and 
we  can  neglect  this  contribution. 
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Considering  the  inhomogeneous  equation,  we  observe  that  we 
could  rewrite  it  as 

-iojt 


mv  + myv  = eE  e 
0 


(29) 


since  u = v and  ii  = v.  This  is  advantageous  since  we  would 

like  to  correlate  j = Nev  to  the  microscopic  process.  A trial 

“iu)t 

solution  V = v^e  gives  us 


and  for  j we  obtain 


D = 


V = 


Ne‘ 


e 

m 


-iu)  + Y 


(30) 


m 


-iu)  + Y 


(31) 


iKx 


Considering  that  j now  also  depends  on  x as  e we  sub- 
stitute for  j from  equation  (31)  and  E = E^e^ into  the 
wave  equation  (26)  , and  obtain 

-K*  = - 


Ne‘ 


-10) 


0 0 


m 


-iu)+Y 


or 


^ 


1 - 


Ne‘ 


Let  us  again  set 


K 


*2 


(n  + ik) 


e m io)-Y  0) 
0 

• 2 + 


(32) 


which,  when  substituted 


o)‘  C‘ 

into  equation  (32)  and  separated  into  real  and  imaginary  parts, 
gives 

Ne*/e  m 

n*  - IT*  - 1 2_  (33) 


0)*  + Y * 


tFor  the  conducting  medium,  the  optical  oarameters  will  be 
written  as  n and  Ic  to  avoid  confusion  with  the  dielectric 
medium  pareuneters. 
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and 


Y 


(34) 


Ne^/e  m 
2 n k = 

0)^  + Y*  (i) 

These  expressions  are  different  from  the  ones  we  previously 
obtained  for  the  dielectric.  For  this  case,  there  is  no  re- 
sonant frequency  involved  and,  consequently,  no  singularity  pre- 
sent when  except  at  (i)=0.  For  we  have  k-+-0  and  n-»l;  which 

means,  for  example  that  x-rays  can  penetrate  metals! 

We  again  use  the  abbreviation  Ne^/e^m  = 0)^^  and  plot  iT  and 
k as  functions  of  to.  From  equations  (33)  and  (34)  we  see  that 
for  Y"*^0  and  to-^O,  n^-]c^  is  negative  and  increasing,  that  is 
|n|  < |kl  and  Ic  has  large  values,  while  for  large  lo  we  have  n=l 

and  k^O.  From  equation  (33)  we  obtain  n=k  for  to=(o  if  y = 0. 

__  _ P 

From  the  equation  for  2nk  we  conclude  that  n is  smaller  than  1. 
Therefore  we  have  qualitatively  for  n and  k,  the  behavior 
sketched  in  Figure  3. 


Fig.  3.  Frequency  dependence  of  refractive  index  n and 
extinction  coefficient  k for  a metal. 
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The  region  where  k becomes  small  and  n approaches  1 is  the 
region  for  which  metals  are  transparent  to  electromagnetic  waves. 
For  large  k,  metals  are  opaque  and  the  metal  reflects  strongly. 
The  wave  can  penetrate  only  a small  distance  into  the  metal  be- 
fore it  is  completely  attenuated.  An  expression  for  the  pene- 
tration depth  is  obtained  by  the  introduction  of  n and  k into  the 
expression  for  the  electromagnetic  wave,  equation  (10),  i.e. 

i [p(n+ik)  x-(i)tl 

E = E e ^ 

0 

“to  . _ M — 

-kpx  i[p;nx-u)t] 

» E^e  e ^ (35) 

Here  2k  ^ is  the  absorption  coefficient  a.  The  depth  x=6  for 
which  the  wave  is  attenuated  to  e“*  of  its  incident  value  is 
called  the  skin  depth  and  is  obtained  from  ^ k6  = 1 or 

6 = — = — (36) 

}cu)  a 


The  optical  behavior  of  the  metal,  in  terms  of  the  above  model, 
depends  upon  the  magnitudes  of  w,  and  y (or  t,  the  relaxa- 
tion time).  For  metals,  typical  values  are  t-IO” ‘ ’s  (y-lO ‘ *s“ * ) 
and  u)p»10*®s“*.  In  the  near  infrared  region  w=10**s”*  so  that 
Y=:uKa)p.  Here  k>n  and  for  large  k,  we  see  that  6 will  be  small. 

Note 

Frequently  in  the  literature  there  is  a different  definition  used 
for  relationships  between  K*  and  n and  k,  that  is 

(K*)  - JUL-  tn(l  + ik)]2 

C* 


- InMl  - IT*)  + 2i  n ITl 

C* 


(37) 
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(38) 


For  the  real  and  imaginary  part  we  would  then  have 


U) 


nMl  - kM  = 1 - 


(jo*  + Y* 


and 


2n  IT  = E- 


u)^  + Y ^ 


(39) 


C.  Absorption  of  Electromagnetic  Waves  by  Molecules 

1.  Normal  Vibrations 

In  considering  the  dynamics  of  a molecule  composed  of 
N atoms,  3N  coordinates  are  required  to  describe  its  trans- 
lational, rotational  and  vibrational  motion.  If  we  exclude 
linear  and  diatomic  molecules  from  our  consideration,  3 
coordinates  are  used  each  for  the  description  of  the  trans- 
latory  and  rotatory  motion  of  the  entire  molecule.  The  re- 
maining 3N-6  coordinates  describe  the  vibrational  motion  of 
the  molecule,  that  is  we  have  3N-6  vibrational  degrees  of 
freedom.  (Linear  and  diatomic  molecules  have  only  3N-5.) 
For  the  determination  of  the  vibrational  frequencies  it  is 
necessary  to  know  the  kinetic  as  well  as  tne  potential 
energy  of  the  atom  in  the  molecular  system.  In  the  Born- 
Oppenheimer  approximation,  one  considers  the  electrons, 
which  are  contributing  to  the  bonding  in  the  molecule,  to 
move  much  faster  than  the  nuclei  (with  or  without  attached 
core  electrons) . The  bonding  electrons  are  responsible  for 
the  potential  energy  of  the  force  field  in  which  the  nuclei 
are  undergoing  their  vibrations. 

The  coordinates  used  for  the  description  of  the  vibra- 
tional problem  may  be  initially  of  various  kinds  such  as 
cartesian  coordinates  or  internal  coordinates.  The  kinetic 
and  potential  energy  is  expressed  in  these  coordinates. 
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The  potential  energy  V is  expressed  as  a Taylor  series 
about  the  equilibrium  position  of  the  molec'  .e 


3N-6 


V = V 


_1_V 
3 c 


(40) 


+ (higher  order  terms) 

where  the  q^'s  are  the  mass  weighted  coordinates  of  the 
various  point  masses  (nuclei  with  and  without  core  elec- 
trons). In  the  harmonic  oscillator  approximation,  terms 
higher  tham  second  order  are  neglected.  Since  0V/Oq^)^ 
is,  by  definition  of  the  equilibrium  position,  equal  to 
zero,  and  is  a scale  factor  (constant) ,the  potential 
energy  becomes  a quadratic  expression  in  i»e. 

3N-6 

2V  « I Jc,  . q.  q.  (40) 

i,j-l  ^ J 

where  the  k^^'s  are  the  force  constants.  Using  these 
coordinates,  the  Icinetic  energy  is  expressed  as 

3N-6 

2T  - I 4 » (42) 

i-1  ^ 


Then,  Lagrange's  formalism  yields  the  equations  of  motion 
3N-6 

qj  + q^  ■ 0 j > 1,  2,  ...,  3N-6  (43) 


A solution  of  this  eigenvalue  problem,  that  is  a solu- 
tion of  the  system  of  linear  second  order  differential 
equations,  is  obtained  by  introducing  a trial  solution  into 
the  above  equations: 


*^i  ” ^i 


(44) 
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Here  u)=2Trf  where  f is  the  vibrational  frequency  to  be  deter- 
mined, and  is  a phase  factor.  The  systems  of  differential 
equations  yield  a set  of  simultaneous,  homogeneous  linear 
algebraic  equations  in  the  3N-6  coordinates. 


3N-6 

I (k..  - 6.  .0)*)  A.  = 0 j = 1 3N-6  (45) 

i = l J 


This  system  of  equations  can  be  solved  only  if  the  secular 
determinant  vanishes,  i.e. 


Sr 


12 

22 


•l,3N-6 


‘3N-6,1 


^3N-6,2  . 


'^3N-6,3N-6"‘*'' 


(46) 


Solution  of  the  polynomial  of  order  3N-6  will  yield  at  most 
3N-6  different  to^'s.  It  may  be  possible  that  some  of  the 
w^'s  are  the  same.  In  such  a case,  the  associated  vibra- 
tional frequencies  are  said  to  be  degenerate.  For  each  of 
the  3N-6  values  of  o)  we  will  obtain  a set  of  A^j^'s  from  the 
system  of  algebraic  equations.  These  sets  of  A^j^'s  are 
called  eigenvectors.  We  can  use  the  3N-6  eigenvectors  to 
compose  a transformation  matrix  L of  the  form 


'll 


l,3N-6 


'3N-6,1 


(47) 


3N-6,3N-6 


/ 
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This  matrix  transforms  the  original  coordinate  system  into 
the  normal  coordinate  system. 


1 

1 1 

^ -1 

= L 

1 : , 

^3N-6  j 

^ ‘^3N-6  J 

The  importance  of  this  transformation  can  best  be  seen  if 
we  go  back  to  the  expression  for  the  kinetic  and  potential 
energy.  Both  are  transformed  by  the  above  matrix  from 
their  dependence  on  the  into  dependence  on  , i.e. 

3N-6  T 3N-6 

2T.  I 

i=l  ^ i-1  ^ 


3N-6 


2V  = 


L JN-6 


(49) 


where  the  k^  are  combinations  of  the  k^^  . The  equation  of 
motion  appears  in  normal  coordinates  as 


Qi  + kf  Qf  - 0 (50) 

that  is,  the  different  equations  are  no  longer  coupled  to 
one  another  and,  therefore,  each  coordinate  is  associated 
solely  with  one  vibrational  mode  of  the  molecule.  In  nor- 
mal modes  all  the  atoms  vibrate  in  phase,  that  is,  they 
all  pass  through  their  equilibrium  positions  at  the  same 
time. 

The  application  of  quantum  mechanics  to  the  system  of 
N atoms  does  not  change  the  results  drastically.  The 
Schroedinger  equation  is  solved  in  terms  of  normal  coordi- 
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nates.  The  modes  of  the  molecular  motion  are  associated 
with  wave  functions  and  eigenvalues.  The  classical  vibra- 
tional frequency  is  the  same  as  the  one  obtained  from 
quantum  mechanics.  However,  selection  rules  are  introduced 
and,  consequently,  the  absorption  intensity  is  modified. 


2 . Selection  Rules  for  a Vibrating  Molecule 

A vibrating  molecule  may  absorb  a light  quantum  hf , 
where  f corresponds  to  the  energy  difference  between  the 
initial  and  exited  (finite)  state.  The  Bohr  frequency  con- 
dition is 


f = 


E,  - E. 
f 1 


(51) 


where  h = PlancJc's  constant.  The  selection  rules  are  ob- 
tained by  expanding  the  electric  dipole  moment  y of  the 
molecule  in  terms  of  normal  coordinates,  i.e. 

3N-6  ( ^ 

y = y + 'I  I Q„  + (higher  order  terms)  (52) 

K=1  ^ 

Again,  in  the  harmonic  oscillator  approximation,  the  higher 
order  terms  are  neglected  for  the  dipole  moment  matrix 
element  / i|^^y  ij^^dT.  We  then  have 


K=1 


'i'fdx, 


(53) 


ij;.  may  be  written  as  a product  of  the  initial  states  of  all 
i th 

vibrations  3N-6  and,  if  only  the  j state  is 

'I'i  = n ip 

singly  exited,  may  be  written  as 


3N-6 

. n ii> 

,j  OK 


(54) 
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Introduction  of  the  wave  function  into  the  expression  for  the 
matrix  element  of  the  dipole  moment  yields 


3N-6 

n I// 

K ^ 


3N-6 

n 

K?<j 


3N-6 


. di 

13 


3N-6 


The  first  integral  vanishes  since  y is  a constant  and 

0 

di  = 0 from  the  orthogonality  condition.  In  the 

second  integral  all  terms  in  the  sum  of  the  type  /\Jj  Q 

OK  K 

are  zero  by  symmetry.  The  only  non-zero  term  is 


This  term  would  be  zero  if 


( 3 S . jo  J*ci  Oj  j 
If  (rao 

1 1 J 


would  vanish  or  if  we  would 
only  one  vibration  to  be 


not  have  initially  chosen  only  one  vibration  to  be 

singly  excited. 

The  general  selection  rule  for  harmonic  oscillator  vibra- 
tional transition  is 

Aj  = ± 1 (57) 

The  permanent  dipole  moment  has  no  influence  on  the  vibrational 
transition,  but  the  change  in  the  dipole  moment  with  respect  to 
the  vibrational  coordinate  (Q^^)  is  important  and  must  be  non- 
zero. The  harmonic  oscillator  approximation  is  sufficient  to 
explain  the  most  intense  transitions.  By  considering  the  in- 
fluence of  the  higher  order  terms  in  equations  (40)  and  (55)  on 
the  dipole  moment  y,  weaker,  previously-forbidden,  transitions 
may  be  predicted. 
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3 . The  Concept  of  Group  Frequencies 

Reviewing  the  near  infrared  absorption  spectra  of  a large 
number  of  chemical  compounds , one  observes  that  absorption 
bands  at  certain  frequencies  are  related  to  the  occurence  of 
certain  molecular  groups  in  the  molecule.  For  example,  the 
bands  around  the  3000  cm” * region  are  related  to  C-H  stretch- 
ing vibrations.  A careful  study  of  this  phenomenon  has  re- 
vealed that  the  relationship  between  observed  bands  and  molec- 
ular groups  is  almost  unique  and  therefore  one  speaks  of  the 
fingerprint  region  if  the  analytical  aspect  of  the  near  infra- 
red region  is  considered.  The  correspondence  between  normal 
vibrations  and  group  frequencies  is,  for  most  of  the  cases, 
simple.  Let  us  consider  the  normal  mode  of  a certain  fre- 
quency which  is  also  listed  as  group  frequency.  The  atoms  be- 
longing to  the  group  frequency  vibrations  are  undergoing  large 
amplitude  vibrations  in  the  normal  mode  compared  to  all  the 
other  atoms.  From  the  point  of  view  of  interaction  with  the 
electromagnetic  freld  the  effect  of  all  the  other  atoms  may  be 
neglected. 

D.  Lattice  Vibrations  in  Dielectric  Media 
1 . Bulk  Material 

The  periodic  arrangement  of  molecules  in  a crystal  lat- 
tice leads  to  new  types  of  absorption  phenomena  which  are 
observed  in  addition  to  the  absorption  due  to  the  internal 
vibrations.  For  most  molecules  in  crystal  form,  lattice  vi- 
bration absorption  is  observed  in  the  middle  and  far-infrared 
regions.  Only  for  small  molecules  may  we  observe  this  type 
of  absorption  in  the  near  infrared.  Consequently,  the  one 
dimensional  linear  chain  may  well  serve  as  a model  for  the 
discussion  of  this  phenomena.  In  this  model,  the  unit  cell 
contains  two  atoms  with  different  effective  charges  so  that 
a contracting  force  exists  between  an  atom  and  its  nearest 
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neighbors.  One  can  see  inunediately  the  occurrence  of  two 
fundamental  modes , one  in  which  the  two  atoms  of  the  unit 
cell  move  in  opposite  directions  and  one  in  which  they 
move  in  the  same  direction.  The  former  is  called  the  longi- 
tudinal optical  mode,  with  frequency  , whereas  the  latter 
is  the  longitudinal  acoustical  mode,  with  zero  frequency.  A 
more  detailed  analysis  leads  to  the  dispersion  relation  of 
the  modes,  that  is,  the  relationship  between  the  frequency 
and  wave  vector  for  all  possible  modes  in  the  Brillouin  zone. 

If  we  extend  these  considerations  to  three  dimensions,  we 
find  that  there  are  two  transvers  optical  and  two  transverse 
acoustical  mode  branches  possible,  in  addition  to  the  longi- 
tudinal one.  Of  all  these  modes,  we  are  only  interested  in 
those  which  lead  to  strong  optical  absorption.  For  inter- 
actions with  an  electromagnetic  field  to  occur,  we  need,  as 
usual,  a component  of  the  electric  field  in  the  direction  of 
the  change  of  the  dipole  moment  of  the  chain.  Since  the 
transverse  modes  are  degenerate,  we  have  only  two  branches 
of  modes  to  consider.  Selection  rules  predict  that  fre- 
quencies belonging  to  zero  wave  vectors  are  the  fundamental 
frequencies.  Difference  amd  summation  frequencies  may  occur 
at  other  non-zero  wave  vectors  with  far  lesser  intensity. 

The  fundamental  frequencies  can  easily  be  visualized  by  con- 
sidering the  chain  as  composed  of  two  sublattices  of  equal 
atoms  having  sublattices  oscillating  rigidly  against  one 
another.  This  results  in  a large  change  in  the  dipole  mo- 
ment amd,  consequently,  in  a large  absorption.  The  absorp- 
tion corresponding  to  this  fundamental  vibration  is  called, 
in  ionic  crystals,  the  reststrahlen  band.  The  names  re- 
late to  the  fact  that  in  very  strong  absorption,  the  incident 
wave  is  attenuated  in  such  a small  surface  layer  that  the 
surface  reflects  most  of  the  radiation,  as  is  found  for  a 
metal . 
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Theoretically  derived  quantitative  expressions  for  the 
absorption  coefficient  are  difficult  to  obtain  since  it  is 
necessary  to  know  the  wave  functions  of  the  ground  and 
excited  state  in  order  to  calculate  the  change  in  the  di- 
pole moment.  Experimentally,  one  proceeds  as  discussed  at 
the  beginning  of  this  section,  and  determines  k,  from 
which  a can  be  obtained. 


2 . Dielectric  Powders 

We  have  just  seen  that  the  fundamental  absorption  pro- 
cess in  ionic  type  crystals  involves  the  transverse  optical 
mode.  In  special  cases,  it  may  also  involve  the  longitudi- 
nal optical  mode.  If  the  diameter  of  the  crystal  becomes 
very  small,  the  frequency  of  the  fundamental  vibration 
shifts  to  higher  values.  It  has  been  observed,  e.g.  in 
MgO  micro-crystals  of  radius  50A  to  800A,  that  the  lattice 
absorption  peak  shifts  up  by  10  cm  The  reason  for  this 
shift  can  be  seen  by  considering  the  number  of  atoms  on  the 
surface  of  the  crystals  (having  neighbors  partially  surround- 
ing them)  compared  to  the  number  of  atoms  in  the  interior  of 
the  crystal  (having  neighbors  completely  surrounding  them) . 

For  a particle  radius  of  50A,  about  25%  of  the  ions  are  on 
the  surface,  whereas  for  a radius  of  380A  there  are  only  3% 
on  the  surface.  Since  the  shift  is  to  higher  frequencies 
with  decreasing  size,  one  may  conclude  that  the  lattice 
spacing  on  the  surface  is  smaller.  This  may  indeed  be  proven. 
As  a consequence  of  this  finding,  one  expects  a strong  in- 
fluence on  the  lattice  vibrational  frequency  of  the  surround- 
ing matrix  due  to  the  polarization  charges.  One  finds,  how- 
ever, a different  frequency  tOp  which  is  related  to  the  trans- 
verse optical  mode  frequency  by 


e + 2e 

2 0 m 2 

F _ . T 


e + 2e 


m 


(58) 
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where  and  are  the  low  and  high  frequency  dielectric 
constants  and  the  dielectric  constant  of  the  surrounding 

medium.  The  frequency  shift  due  to  this  effect  can  be  re- 

« 

markable.  In  BaTiO.  , for  example,  the  lattice  vibrational 

^ -1 

frequency  of  the  bulk  material  is  at  33  cm  , whereas  for 

small  spheres  in  powder  form  in  a surrounding  medium  with 

e “1,  one  observes  a value  of  400  cm 
m 

Experiments  or.  small  particles  are  hampered  by  the  scat- 
tering effects  of  the  particles.  Quantitative  experimental 
values  can  be  obtained  via  the  dielectric  constant  and  from 
that,  the  determination  of  a may  be  made. 

E.  Absorption  of  Electromagnetic  Waves  by  Semiconductors 

For  absorption  by  semiconductors  in  the  near  infrared 
spectral  region  we  consider:  (a)  Transitions  of  electrons 
across  the  band  gap,  (b)  free  carrier  absorption  and  inter- 
band absorption,  and  (c)  multiphonon  absorption. 

1.  Band  Gap  Absorption 

The  strongest  absorption  in  the  near  infrared  is  due  to 
the  transition  of  electrons  from  the  valence  band  to  the  con- 
duction band.  These  tramsitions  are  only  possible  if  the 
energy  of  the  absorbed  photon  is  comparedjle  or  larger  than 
the  band  gap  energy,  that  is  hf  i E^  . For  photons  with  hf 
smaller  than  E^  these  transitions  are  not  possible.  A plot 
of  the  absorption  coefficient  a(f)  against  the  photon  energy 
shows  therefore  a strong  increase  of  absorption  near  hf=Eg; 
this  is  called  the  band  edge  absorption.  For  hf  < E^  the 
material  is  relatively  transparent.  The  absorption  in  the 
hf  > Eg  region  is  relatively  strong;  a has  values  of  the 
order  of  10 "/cm. 

The  transition  from  the  valence  to  the  conduction  band 
can  be  a "direct"  or  an"indirect "transition , depending  on 
the  structure  of  the  bands.  The  energy  difference  E^ 
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between  the  bands  depends  in  general  on  the  crystal  struc- 
ture. Direct  transitions  occur  if  the  minimum  of  the  con- 
duction band  and  the  maximum  of  the  valence  band  both  are 
situated  at  the  same  K value  (e.g.  at  K = 0).  The  electron 
then  makes  a "vertical  transition"  without  change  in  mo- 
mentum (K  is  conserved)  and  the  energy  difference  is  = hf. 
An  indirect  transition  occurs  if  the  minimum  of  the  conduc- 
tion band  and  the  maximum  of  the  valence  band  are  not  situ- 
ated at  the  same  point  in  K space.  In  order  to  conserve  mo- 
mentum a phonon  is  either  absorbed  or  emitted.  Since  the 
photon  momentum  is  small  in  the  energy  region  under  consid- 
eration, the  phonon  makes  up  for  the  change  in  K.  On  the 
other  hand,  the  phonon  energy  is  small  compared  to  the  energy 
of  the  absorbed  photon  so  that  hf  = E^.  The  absorption  pro- 
bability for  an  indirect  process  is  smaller  that  that  for  a 
direct  process. 

(a)  Absorption  coefficient  for  direct  transitions 

Let  us  assume  a parabolic  relationship  between  E and 
K at  the  maximum  of  the  valence  band  and  at  the  minimum 
of  the  conduction  band  where  the  direct  transition  occurs. 
The  absorption  coefficient  is  then  found  to  be: 


a (f)  = A / (hf  - E ) 

g 


where 


3/2 


n C h*  m_* 


(59) 


(60) 


q is  the  charge  of  the  electron,  and  m^*  are  the 

effective  masses  of  holes  and  electrons  respectively,  n 
is  the  real  refractive  index,  C is  the  speed  of  light, 
and  h is  Planck's  constant.  For  n - 4 and  mj^*  = m^*  = m^  , 
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(61) 


where  is  the  free  electron  rest  mass,  we  have 
a(f)  - 2 X 10“  (hf  - E 

g 

a(f)  is  given  in  cm~ * if  h and  are  inserted  in  units  of 
electron  volts. 


(b)  Absorption  coefficient  for  forbidden  direct 
transitions 

Sometimes  direct  transitions  are  forbidden  by  selec- 
tion rules  for  K = 0 (symmetry) , but  allowed  for  K ^ 0. 
For  this  case  the  absorption  coefficient  is  found  to  be: 


where 


a(f)  = A'  (hf  - Eg)’/* 


''‘'h*”*e* 

q*  -f  m^* 

n C h*  m^*  nij^*  (hf ) 


(62) 


(63) 


and  for  the  case  discussed  above  [paragraph  (a) ] we 
have 

^ (hf  - E^)’/* 
alf)  = 1.3  X 10“ a 


hf 


cm 


- 1 


(64) 


again. 


hf  and  E are  in  ev's. 

g 


(c)  Absorption  coefficient  for  indirect  transitions 

The  indirect  transitions  are  two-step  processes  and, 
consequently,  the  interaction  of  photons  with  electrons 
and  phonons  have  to  be  considered.  The  absorption  co- 
efficient has  been  shown  to  be 


a(f) 


(hf  - E - E^) 

— a E_ 


1+  exp(-Ep/)cT) 


-] 


(hf  - . Ep) 

exp(-E  AT) 


9 


(65) 
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where  is  the  energy  of  the  phonon.  The  first  term  is 
due  to  a process  with  phonon  absorption,  the  second  with 
phonon  emission.  The  constant  A has  a complicated  ap- 
pearance. In  many  cases  it  is  determined  by  an  experi- 
mental  fit  to  the  data. 


2 • Free  Carrier  Absorption 

Free  carrier  absorption  is  of  importance  in  semicon- 
ductors as  it  is  in  conductors.  Absorption  due  to  free  car- 
riers may  be  observed  in  the  near  and  far  infrared.  The  ab- 
sorption of  photons  by  the  electrons  in  the  periodic  field 
of  the  crystal  is  strictly  not  allowed  as  a consequence  of 
conservation  of  momentum.  However,  a perturbation  involving 
lattice  vibrations  (phonons)  relaxes  the  restrictions  and  ab- 
sorption of  photons  is  possible.  The  absorption  coefficient 
may  be  calculated  on  the  basis  of  classical  theory  by  assum- 
ing an  oscillating  electron  without  a restoring  force  and 
using  an  effective  mass  m^*  in  place  of  the  electron  rest 

mass  m . One  obtains: 
e 


n 


- - e = 


w*  + Y* 


(66) 


NqVm  *e 

2nka)  = - — — (67) 

u)*  + 


where  y = — ^ and  t is  the  collision  time.  For  wave- 
lengths up  to  a few  hundred  micrometers , the  absorption  co- 
efficient may  be  expressed  as: 


(68) 
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where  \ is  the  wavelength;  q is  the  electron  charge;  is 
the  electron  mass;  is  the  electron  mobility;  C is  the 
speed  of  light;  n is  the  refractive  index  and  is  the 
static  dielectric  constant.  The  order  of  magnitude  of  the 
absorption  coefficient  a(f)  is  10  - 100  cm“ ‘ for  free  car- 
rier ad5sorption  in,  for  example,  n-type  InAs  in  the  3-15ym 
region.  This  value  is  two  orders  of  magnitude  smaller  than 
that  for  band  edge  edjsorption.  Interband  edjsorption  and  ab- 
sorption due  to  holes  results  in  absorption  coefficients  of 
the  same  magnitude. 

3 . Multiple  Phonon  Absorption 

The  absorption  due  to  lattice  vibrations  occurs,  for 
most  material,  in  the  middle  and  far  infrared.  For  ionic 
crystals  these  lattice  vibrations  show  up  as  reststrahlen 
bands.  These  bemds  are  due  to  fundamental  modes  of  the  lat- 
tice with  K * 0,  K being  the  crystal  momentum  vector.  Mul- 
tiple phonon  absorption  of  the  summation  type  can  be  ob- 
served in  the  near  infrared  region.  Absorption  due  to  such 
processes  are  relatively  weaJc,  because  these  processes 
second  order  effects.  They  appear  as  modulation  on  top  of 
the  free  carrier  absorption. 

F.  Absorption  of  Electromagnetic  Waves  by  Metals 

1.  Bulk  Materials 

The  dependence  of  the  optical  constants  on  the  micro- 
scopic parameters  of  bulk  conducting  materials  has  been  dis- 
cussed in  the  first  chapter  on  "simplified  dispersion  theory". 
In  this  simple  approach  the  plasma  frequency  and  relaxation 
time  are  the  characteristic  parameters  for  the  atbsorption 
process.  This  process  is  similar  to  the  2Ubsorption  in  semi- 
conductors as  we  have  just  seen  in  section  I.  E.2  "free 
carrier  absorption". 
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2.  Metal  Powders 


By  considering  metal  powder  instead  of  the  bulk  mat- 
erial one  anticipates  changes  analogous  to  those  observed 
for  dielectric  bulk  material  and  their  powders.  The  impor- 
tant parameter  is  the  plasma  frequency  and  for  particles  of 
about  lOA  diameter,  a broadening  of  the  plasma  resonance 
peak  is  observed.  Theoretically,  this  broadening  is  des- 
cribed by  a decreased  electron  mean  free  path  due  to  scat- 
tering at  the  particle  surface.  This  macroscopic  explana- 
tion has  been  more  successful  than  considerations  from  a 
quantum  mechamical  point  of  view.  In  thin  conducting  films 
on  the  other  hand,  one  has  observed  effects  which  are  inter 
preted  as  the  breakdown  of  the  continuous  conduction  band 
into  discreate  states  by  reducing  the  thickness  of  the  film 
from  bulk  material  dimensions.  The  evaluation  of  the  ab- 
sorption coefficient  from  experimental  measurements  pro- 
ceeds as  outlined  before. 

G.  Absorption  of  Electromagnetic  Waves  by  One-Dimensional 
Conductors 

Recently  a number  of  linear  conducting  materials  have 
been  developed.  These  polymeric  materials  have  metallic  be- 
havior in  the  direction  of  the  chain  and  almost  insulator-like 
properties  in  the  direction  perpendicular  to  the  chain.  An  ex- 
ample of  a polymeric  one-dimensional  conductor  (not  recognized 
as  such  initially)  is  the  H-sheet  polarizer  invented  by  E.  Land 
in  1938.  This  polarizer  is  made  of  polyvinyl  alcohol  chains 
which  are  aligned  by  uniaxial  stretching  of  the  material.  The 
sheet  is  impregnated  with  Iodine  by  immersing  it  into  a special 
ink.  The  iodine  atoms  or  molecular  clusters  are  believed  to  do- 
nate electrons.  These  electrons  constitute  the  free  charge  car- 
riers involved  in  the  conduction  process  as  in  a metallic  con- 
ductor. Incident  light  with  its  electric  field  vector  oriented 
parallel  to  the  polymer  chain  is  absorbed  so  that  the  sheet  acts 
as  a polarizer. 
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As  a result  of  the  considerable  interest  in  linear  conductors 
a number  of  one-dimensional  linear  polymeric  conductors  have  been 
extensively  studied,  e. g.  (SN)^  and  polyacetylene-iodide.  These 
materials  are  black  and  light-weight,  emd  for  polyacetylene-iodide 
the  absorption  in  the  near  infrared  region  is  broad  and  strong, 
but  without  distinctive  eJssorption  features.  A quantitative  eUa- 
sorption  study  is  underway  in  our  laboratory. 

II.  QUANTITATIVE  ANALYSIS  OF  ABSORBING  MATERIALS 

An  electromagnetic  wave  incident  upon  a material  will  inter- 
act with  it,  and  have  its  energy  changed,  by  three  types  of  mech- 
anisms; (1)  Reflection,  in  which  a fraction  of  the  incident  radi- 
ation is  deflected  backwards;  (2)  scattering,  in  which  a fraction 
of  the  incident  radiation  is  scattered  in  all  directions  (charac- 
teristics depend  upon  relationship  between  particle  size  and  wave- 
length of  radiation,  optical  inhomogeneities,  etc.);  and  (3)  ab- 
sorption, the  general  term  describing  the  conversion  of  radiant 
energy  into  thermal  energy.  Generally,  all  these  processes  are 
present  to  varying  degrees  aind,  are  sensitive  to  the  wavelength 
of  the  incident  radiation.  We  are  interested  here  in  the  ab- 
sorption process  and  assume  that  the  beam  incident  normally  upon, 
and  entering  the  medium,  has  a radiant  power  I^.  Upon  passing 
through  a uniform  layer  b units  thick,  the  radiant  power  is 
attenuated  to  a value  I (b)  given  by  the  absorption  (Lambert- 
Bouguer)  law 

Kb)  - I e (69) 

0 

The  product  a(v)b  Is  the  optical  density  or  absorbance  and  a(v) 
is  the  (frequency-dependent;  “ y)  absorption  coefficient.  As 
mentioned  in  paragraph  I.  A,  a(v)  is  related  to  the  Imaginary 
part  of  the  complex  index  of  refraction  k(v)  » ImIn(v)*I  by 

a(v)  « —1^  k(v)  ■ 4iTvk(v)  (70) 
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The  extinction  coefficient  k(v)  is  a function  of  frequency  [as 
is  n(v)l.  Equations  (69)  and  (70)  show  that,  in  principle,  one 
has  only  to  measure  the  total  Scunple  thic)cness  b and  to 

extract  values  of  a(u)  and,  if  desired  )c(v),  for  material^  in 
the  form  of  uniform  slabs.  If  the  material  of  interest  is  dis- 
solved in  a suitable  (i.e.,  non-absorbing,  non-interacting,  etc.) 
solvent,  Eqn.  (69)is  written  as 

Kb)  = I^e  • (71) 

usually  referred  to  as  Beer's  Law.  Here  b is  the  optical  path 
length  through  the  solution  of  concentration  c moles/liter 
characterized  by  a molar  extinction  coefficient  £ (v)  measured 
in  units  of  liter/mole-cm.  Clearly  for  a liquid  film  confined  to 
a cell  of  thickness  b cm,  comparison  of  equations  (69)^(70)  and 
(71)  gives 

e (V)  = vk(v)  (72) 

The  concentration  may  be  obtained  from 

c = -g-  d.  (73) 

where  p is  the  solute  density  (in  gm/m£),  M is  the  molecular 
mass  (in  gm/mole)  , and  di  is  volume  fraction  of  solute  [i.e.,  m£. 
(solute) /liter  (solution)].  Clearly,  in  assuming  Beer's  Law 
to  hold  for  a solution,  we  must  determine  experimentally  the 
absorbance  A and  either  e(v)  or  c.  In  routine  analysis  of  mat- 
erials, e (v)  is  usually  known  and  the  concentration  c is  the 
desired  quantity. 

A quantity  of  considerable  interest,  especially  as  related 
to  the  study  of  molecular  dynamics,  is  the  integrated  absorption 
intensity,  defined  here  as 

/V2 

a(v)dv  (74) 

Vl 
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Here  v and  v are  frequencies  (wave  numbers)  sufficiently  far 
1 2 

removed  from  the  band  center  that  a(v  ) = a(v  ) =0. 

1 2 

Theories  relating  the  microscopic  processes  to  a(v)  are 
mathematically  quite  involved.  The  resonance  absorption  of  a 
molecule  in  the  gaseous  state,  for  example,  depends  on  both  the 
vibrational  frequency,  and  the  frequency  of  the  overall  rotation 
of  the  molecule.  Interactions  between  the  vibrational  and  rota- 
tional motion  occur.  The  analysis,  theoretical  or  empirical, 
should  reveal  the  observed  frequencies.  The  intensities  may  then 
be  calculated  by  assuming,  for  example,  a Boltzmann  equilibrium 
distribution  of  the  oscillators.  The  effect  of  the  collisions 
between  the  molecules,  that  is,  line  broadening,  might  be  con- 
sidered as  well  as  other  factors  affecting  the  microscopic  pro- 
cesses, and  therefore  the  absorption  coefficient.  Calculations 
are  only  possible  for  the  simplest  of  molecules.  For  larger 
molecules  one  needs  interaction  constants  of  the  vibration- 
rotation  interaction  or  the  vibrational  wave  functions  of  the 
ground  and  first  excited  state.  These  are  difficult  to  calculate 
and  enpirically  not  easy  to  determine.  As  an  "over-simplified" 
rule,  however,  in  most  cases  the  absorption  coefficient  is  pro- 
portional to  the  resonance  frequency  and  the  square  of  the  matrix 
element  of  the  dipole  moment.  In  liquids,  the  overall  rotation 
of  the  molecule  is  suppressed  and  the  collisions  between  molecules 
play  a more  dominant  role.  The  edjsorption  process  might  be  con- 
sidered as  consisting  of  resonance  eibsorption  and  orientational 
ad^sorption.  The  former  is  an  absorption  centered  about  the  vi- 
brational resonance  frequency  where  the  rotational  states  do  not 
exist  as  individual  levels,  but  contribute  to  a broad  absorption 
band.  The  latter  depends  on  the  orientation  (and  reorientation) 
of  permanent  dipoles  in  an  applied  electromagnetic  field  (Debye 
absorption) . An  empirical  description  of  these  two  processes 
is  of  the  form 

o(t)  « a e cos  (bit  + <|»)  + a e 

1 2 
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If  the  constants  a , t,  a and  (}>  are  determined,  a(v)  may  be  ob- 

1 2 

tained  through  a Fourier  transformation.  Here  t is  a character- 
istic relaxation  time  associated  with  the  orientation.  Theoret- 
ical calculations  of  a(v)  using  procedures  from  statistical  mech- 
anics, as  well  as  quantum  mechanics,  are  difficult  euid  lengthy. 
However,  this  is  a very  active  field  since  one  can  deduce  the 
correlation  function  and  obtain  em  estimate  of  the  molecular  be- 
havior subsequent  to  an  assumed  initial  configuration.  The  time 
interval  considered  is  of  the  order  of  10“*‘-10”‘^  sec. 

Experimental  Methods 

The  use  of  the  Lambert-Bouguer-Beer  Law  appears  deceptively 
straightforward.  While  the  method  is  not  sophisticated,  there 
are  certain  precautions  that  must  be  considered  so  that  the  data 
yield  true  material  parameters  rather  than  experimental  artifacts. 
We  will  consider  in  limited  detail,  the  more  significant  of  these 
features . 

A.  Liquids 

1.  Deviations  from  Beer's  Law  and  Photometric  Error 

The  attenuation  of  a beam  of  light  as  it  passes  through 
a specimen  depends  upon  the  samples  thickness  b,  its  ab- 
sorption coefficient  a and  the  particular  wavelength  X or 
frequency  (actually  wavenumber  v = X~ ' ) . Ultimately,  the 
attenuation  depends  upon  the  number  of  absorbing  particles 
(molecules)  "seen"  by  the  optical  beam.  Provided  that  the 
number  of  particles  and  the  volume  containing  them  remain 
constant,  the  absorption  coefficient  a(v)  suffices  to  char- 
acterize the  substance.  Such  is  the  case  for  solids  and 
pure  liquids.  For  solutions  or  mixtures  of  liquids,  however, 
this  condition  is  not  fulfilled  and  it  becomes  necessary  to 
relate  the  adasorption  to  the  sample  concentration.  (The 
analogous  situation  for  gases  would  require  that  the  effect 
of  pressure  be  considered.)  ' iiis  is  done  in  the  case  of 
Beer's  Law,  since  the  absorbance  is  given  by 
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A = e (v)  c b (75) 

For  a substance  to  "obey"  Beer's  Law  then,  the  absorbance 
must  be  proportional  to  the  concentration.  Unfortunately, 
this  is  not  always  found  and  deviations,  both  real,  as  in 
the  case  of  refractive  index  changes  due  to  high  solute  con- 
centrations, and  apparent,  due  to  instrumental  limitations 
and  non-symmetrical  chemical  equilibrium,  are  encountered. 

The  real  deviations  arise  when  the  measured  absorption  is  no 
longer  the  sum  of  the  absorptions  for  all  separate  particles. 
More  specifically  the  relation 

N N 

A >2  Ai  - b 
i=l 

where  i labels  the  species  and  N is  the  total  number  of 
species  present,  is  not  obeyed.  This  non-ideal  behavior  is 
a result  of  formation  of  molecular  aggregates  or  clusters  by 
solute-solute  and  solute-solvent  interactions.  This  type  of 
deviation  may  be  used  to  study  molecular  interactions,  and  is 
therefore,  not  without  some  value. 

The  experimental  deviations  arise  because  Beer's  Law  is 
rigorously  correct  only  for  monochromatic  radiation.  The 
finite  slit  width  of  the  instrument  monochromator  allows  a 
band  of  frequencies  to  be  detected  so  that  the  magnitude  of 
the  deviation  defends  on  the  ratio  of  the  slit  width  to  the 
absorption  banu  'dth  observed.  This  effect  has  no  compen- 
sating benefi*^s,  but  must  be  tolerated.  However,  empirical 
determination  of  actors  affecting  the  absorbance  may  still 
be  evaluated  provided  the  aforementioned  instrumental  fac- 
tors are  maintained  constant. 

The  above  considerations  of  deviations  leads  to  some 
practical  rules-of-thumb  for  choosing  the  experimental  para- 
meters. To  eliminate  clustering  effects  in  solutions,  con- 
centrations below  about  .OlM  should  be  used.  However,  the 
instruments  €dDility  to  detect  small  changes  in  radiant  power 


I 


(v)c. 


(•'6) 
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limits  the  dilution  to  about  10~’m.  Conconunitant  with  the 
desire  for  low  concentrations  is  the  instrument's  need  for 
sufficient  energy  at  the  detector  so  that  photometric  error 
may  be  minimized.  Here  again  a suitable  balance  must  be 
sought  since  for  a low  solute  concentration,  a weak  ab- 
sorption band  results  and  the  instrumeat  "sees"  a small 
difference  between  incident  and  transmitted  power,  and  thus 
a relatively  large  error  in  absorbance.  Alternatively,  for 
a more  concentrated  solution,  the  absorption  band  is  strong 
and  insufficient  power  reaches  the  detector,  also  giving  in- 
accurate results.  An  appropriate  intermediate  range  of  ab- 
sorbance can  be  obtained  starting  with  Eqn.  (75)  rewritten 
as 

* -BS  - -E  ^ - -fe 

Here  we  have  used  the  identity  A = -log  T,  where  T is  trans- 
mission, usually  associated  with  1/1^-  Our  aim  is  to  mini- 
mize the  error  in  the  molar  extinction  coefficient  e(v)  by 
minimizing  the  error  in  T.  The  fractional  error  in  e (v)  is 


de  (v) 

~TNT 


log  e 
e (v ) be 


dT 


(78) 


where  log  means  log  ("<ln"  is  reserved  for  log  ) log  e = 
1/2.30  = .434,  and  dT  is  the  uncertainty  in  T.  To  minimize 
dt/e , differentiate  Eqn.  (78)  and  equate  the  result  to  zero. 
The  result  gives  log  T + log  e = 0,  log  T = -.434,  or  T^^^  = 
.368  for  the  optimum  value  of  transmission  to  minimize  in- 
strument photometric  error. 

To  summarize,  optimum  accuracy  in  the  determination  of 
e(v)  is  obtained  when  the  solute  concentration  is  such  that 
(1)  molecular  clustering  is  avoided  and  (2)  the  transmission 
is  in  the  range  of  30-40%  for  absorption  band  minima. 
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2.  Cell  Losses 


In  stating  the  Lambert-Bouguer-Beer  Law  [Bqns.  (69)  and 
(71)1  the  quantity  is  taken  as  the  incident  beam  power  and 
I the  transmitted  power.  Since  a beam  of  light  is  partially 
reflected  at  the  interfaces  between  two  dielectric  media, 
the  incident  power  arriving  at  the  first  surface  of  a liquid 
cell  is  not  totally  transmitted  to  the  second  surface,  where 
additional  reflection  also  takes  place.  In  fact,  multiple 
reflections  will  take  place  at  all  four  interfaces  of  a liquid 
cell  consisting  of  two  (assumed)  flat  windows  separated  by 
an  (assumed)  uniform  space.  Ignoring  the  convergence  of  the 
optical  beam  in  the  sample  compartment  of  most  infrared 
spectrometers,  the  fraction  of  incident  radiation  reflected 
at  normal  incidence  can  be  calculated  from 


I 

T 


refl 

0 


n 

- n 

i 

n 

(79) 


For  an  NaCl  window  (n^  = 1.54)  in  air  (nj  = 1.00), 

I -ri/I  = .045  or  4.5%  of  the  incident  power  is  reflected 
from  the  first  surface.  Clearly,  for  accurate  quantitative 
work,  the  various  cell  losses  must  be  accounted  for.  One 
method  suggested  by  Draegert,  et  al  [J.O.S.A.  64,(1966)] 

and  used  in  this  reported  work  is  to  consider  the  actual  in- 
cident beam  to  be  represented  by  and  the  transmitted 

beam  by  rl.  The  coefficients  r and  r^ , accounting  for  the 
cell  losses,  need  not  be  equal. 

The  measured  transmission  of  the  cell-solution  system  is 

T = rl/r  I while  the  true  transmission  is  T » e = 

rc  0 0 

I/I^  . Equation  (71)  now  becomes 

UnT^  = An  (-~)  - e(v)cb  (80) 

m r _ 
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I 

If  one  makes  a series  of  measurements  of  T for  various  cell 

m 

thicknesses  b,  holdinq  the  concentration  c constant  and 
plots  in  vs  b,  the  resulting  curve  is  a straight  line 
whose  slope  is  -e(v)c.  The  molar  extinction  coefficient 
£.(v)  can  then  be  extracted  from  the  data  by,  viz.  a least 
squares  fit.  The  ordinate  intercept  of  this  curve  gives 
r/r^ , but  is  of  little  practical  significance.  As  indicated 
above,  we  have  used  this  method  to  obtain  e(v)  for  tri-n- 
butyl  phosphate.  As  a check  on  the  validity  of  the  method, 
we  determined  the  molar  extinction  coefficient  e (v)  and  from 
it,  the  extinction  coefficient  k(v)  = (c/4Tr)  e(v)/v  [see 
Eqn.  (72)]  for  pure,  spectral  quality  chloroform  (CHCl^) • 

Our  results  at  3020  cm” M3 . 31y ) [k (3020)  = 0.011]  agree  very 
well  with  those  reported  by  Hauranek  and  Jones  [Spectrochim. 

^ Acta  32A,  111  (1976)]  [k(3020)  = 0.0115]  using  a considerably 

r 

more  involved  procedure. 

3.  Other  Corrections 

Using  detailed  computer  calculations  for  measured  and 
simulated  data,  Jones  and  coworkers*  have  considered  fac- 
tors contributing  to  spectrophotometric  qrrors  in  determin- 
ing the  optical  constants  of  materials  confined  to  thin 
cells  (l-lOOp).  For  information  only,  we  list  these  factors 
with  our  parenthetical  comments  where  appropriate. 

(a)  Construction  of  precision  (thin)  cells. 

(b)  Distortion  caused  by  interference  and  reflection 
within  the  cell  and  at  cell  surfaces. 

(c)  Polarization  discrimination  of  monochromator. 

(d)  Distortion  caused  by  wedge-shaped  cell,  i.e.,  non- 
parallel windows. 


* 

See  Spectrochima  Acta  32A,  75-123  (1976) 
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(e)  Non-normal  beam  incidence  due  to  beam  convergence  - 
(effects  interference  fringe  pattern  when  b>100ij)  . 

(f)  Surface  roughness  of  window  faces  - negligible  in 
mid  and  far-infrared,  i.e.  for  A - 2.5p  (cell  win- 
dows should  be  demounted  and  repolished  ~ every 
three  months) . 

(g)  Indeterminacy  in  cell  thickness  value  caused  by 
cold  flowing  of  window  material  and  (e)  (cell 
spacing  should  be  measured  before  and  after 
spectrum  is  obtained) . 

(h)  Dispersion  distortion  of  aJasorbing  solutes  in  non- 
absorbing solvents  [important  for  high  concentra- 
tions (i.e.  c ~ 50  mol%)  in  thin  cells  (i.e. 

b ~ ly)]. 

(i)  Various  computational  errors  generated,  for  ex- 
ample, by  truncation  of  wave  number  interval  in 
computer  program  (not  a factor  in  this  work) . 

Many  of  the  above-listed  factors  are  of  significant 
importance  only  if  extremely  accurate  (i.e.  archival)  de- 
terminations of  optical  constants  are  sought.  In  these 
cases  considerable  experimental  control  and  computer  time 
are  required. 

B.  Solid  Films 

1.  Reflection  Losses 

Since  free  standing  solid  films  have  only  two  inter- 
faces, as  compared  to  four  in  a liquid  film  cell,  the  opti- 
cal problems  involved  are  somewhat  simpler.  Reflections 
will  still  occur  at  the  first  and  second  surfaces  and  inter- 
ference fringes  may  appear  if  the  film  thickness  has  the 
appropriate  value.  Fringes  may  be  eliminated  by  (a)  casting 
a film  directly  on  a transparent  window  or  (b)  coating  the 
film  with  Nujol  (mineral  oil)  and  pressing  it  against  a 
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thick  window  [Latinski,  C. , Anal.  Chem.  2071  (1958)]. 

Of  course  if  the  film  is  thick  enough,  the  fringe  separa- 
tion is  less  than  the  monochromator  resolution,  and  inter- 
ference fringes  are  not  observed. 

Reflection  losses  at  the  air/film  interfaces  may  be 
treated  in  a manner  similar  to  that  described  above  for 
liquid  film  measurements.  That  is,  the  spectrum  is  scanned 
repeatedly  using  films  of  varying  thickness.  One  then  uses 
the  equation 


tn  T = 
m 


a (v)  b 


(81) 


to  extract  the  absorption  coefficient  a(v).  The  extinction 
coefficient  k(v)  = Im[n*(v)]  may  then  be  obtained  (see  Eqn. 

(70) ] , if  so  desired. 

2.  Other  Errors 

The  discussion  presented  adsove  in  paragraphs  A. 2 and 
A. 3 are  generally  applicable  for  thin,  free-standing  films, 
with  some  minor  modifications  possibly  required.  Therefore, 
they  will  not  be  further  addressed  here.  Of  greater  impor- 
tance here  are  two  experimental  difficulties  related  to  sample 
preparation.  First;  Is  the  (solid)  sample  available  in  thin 
film  form?  Here  "thin"  means  sufficiently  transparent  in  the 
wavelength  region  of  interest  so  that  photometric  errors 
(see  paragraph  A.l)  are  minimized,  i.e.  so  that  T = 30-40%. 
Second:  Is  the  thin  film  of  uniform  thickness? 

If  the  solid  material  is  insoluble,  such  as  with  some  poly- 
mers, films  may  be  produced  by  rolling  or  pressing  the  bulk 
sample  in  a molten  state.  Both  methods  offer  positive  and 
negative  features.  For  example,  hot  pressing  of  polymer 
beads  on  powder  may  not  offer  sufficient  temperature  control. 


allowing  for  thermal  degradation,  or  may  allow  air  pockets 
or  bubbles  to  form  within  the  film,  causing  photometric 
inaccuracies.  Rolling  samples  may  avoid  these  problems, 
but  introduces  orientation  effects,  possibly  undesirable. 
Sectioning  thin  slices  is  possible  if  a microtome  is  avail- 
able and  if  the  material  is  rigid  enough. 

For  soluble  materials,  casting  films  from  solutions  is  a 
practical  method.  Numerous  techniques  have  been  developed 
to  produce  uniformly  thick  films,  mamy  of  which  have  been 
reviewed  by  Mano  and  Durao  [J.  Chem.  Ed.,  228  (1973)]. 

The  choice  of  solvent,  environment  of  cast  solution  and 
condition  of  the  substrate  surface  are  among  the  relevant 
factors  to  be  considered  here. 


f III.  EXPERIMENTAL  RESULTS 

A.  Polymethylmethacrylate  (PMMA) 

PMMA  samples  were  obtained  from  cast  sheets  of  Lucite 
supplied  by  the  Plastic  Products  and  Resins  Department  of 
E.  I.  DuPont  de  Nemours  6 Co.  Lucite  pieces  were  dissolved 
in  reagent  grade  toluene.  The  resultant  solution  was  poured 
into  teflon  casting  dishes  which  were  set  upon  a leveled 
microflat.  The  rate  of  solvent  evaporation  was  controlled 
by  enclosing  the  casting  dishes  within  a large  glass  con- 
tainer (inverted  fish  temk) . Sample  thickness  was  controlled 
by  varying  the  volume  of  solution  into  the  casting  dish. 

In  spite  of  the  edsove  precautions  undertaken  to  obtain 
uniform  films,  the  resulting  sample  thicknesses  exhibited 
considerable  variation.  The  three  PMMA  films  used  to  deter- 
mine a(v)  were  measured  and  found  to  have  the  following 
values . 
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0.375 

± .029 

mil  = 

9.53 

0.74m 

( 8% 

0.225 

± .050 

mil  = 

5.72 

+ 

1.27m 

(22%) 

0.169 

+ .067 

mil  = 

4.29 

+ 

1.70m 

(40%) 

Film  thickness  was  measured  with  a Starrett  Gauge  (No.  25- 
209) . The  average  of  b "seen  by"  the  optical  beam  is 
listed  above  together  with  the  (%)  deviation.  This  magni- 
tude of  the  deviation  in  film  thickness  is  the  principal 
source  of  error  in  the  determination  of  a(v)  for  PMMA  and  is 
a consequence  of  the  need  to  have  films  thin  enough  so  that 
sufficient  radiant  energy  is  transmitted  and  reaches  the  de- 
tector of  the  spectrophotometer  (see  paragraph  II  ) . Since 
strongly  absorbing  materials  are  of  primary  interest  in  this 
project,  the  problems  associated  with  producing  "uniform  and 
thin"  films  cannot  be  avoided. 

The  infrared  spectra  of  the  three  films  of  PMMA  were  re- 
corded from  4000  to  800  cm" * (2.5  - 32. 5u).  The  value  of 

a(v)  was  calculated  from  Eqn.  (81  ) at  20  cm"'  intervals 
from  4000  to  2000  cm” * (2.5  - 5u) , and  at  10  cm"*  intervals 
from  2000  to  850  cm”*  (5.50  - 12.5vi).  All  values  reported 
were  calculated  by  a linear  regression  analysis  on  the  data 
points  (2.n  T vs  b)  . The  results  are  presented  in  Fig.  4 a 
and  4 b. 

B.  Tri-n-Butyl  Phosphate 

Solutions  of  tri-n-butyl  phosphate  (TBP)  were  prepared  by 
mixing  the  commercially  available  solute  (Fisher  No.  B-404) 
with  spectral  quality  solvents  (CCl^  and  CS^).  Varying  con- 
centrations were  produced  by  measuring  a known  quantity  of 
TBP,  using  Ostwald-Van  Slyke,  Class  A-type  pipets,  into  cali- 
brated volumetric  flasks  and  filling  to  level  with  solvent. 
Precalibrated  stainless  steel,  sealed  liquid  cells  (Inter- 
national Crystal  Laboratories  Model  SL-3)  with  Nad  windows 
and  cell  spacings  from  .025  to  .500  mm  were  used  to  investi- 
gate the  solution  spectra.  No  compensating  solvent-containing 
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cell  was  used  in  the  reference  beam,  as  recommended  by 
Hawranek,  et  al  iSpectrochimia  Acta  32A,  75  (1976)].  Spectral 
quality  CCl^  was  used  as  the  solvent  for  the  4000  - 1300  cm"' 
region  and  CSj  was  used  below  1300  cm”'.  Spectra  were  ob- 
tained with  instrument  parameters  chosen  to  optimize  sensi- 
tivity. 

As  discussed  above  in  paragraph  II  , in  the  Bouguer-Beer 
Law  (Eqn.71  ),  the  value  of  the  molar  extinction  coeffici- 

ent e (v)  should  be  independent  of  concentration.  Values  of 
c < 10” are  recommended  for  solution  measurements,  the  low 
concentration  limit  being  determined  by  instrumental  condi- 
tions. However,  for  real  systems,  one  expects  deviations 
from  ideal  Bouguer-Beer  behavior  due  to  refractive  index 
changes  associated  with  concentration  effects.  Allowing  for 
this  possibility,  solutions  with  concentrations  of  0.0073  M, 
0.0146  M,  0.0292  M and  0.073  M were  prepared.  For  each  con- 
centration, spectral  absorption  measurements  were  made  at 
three  optical  path  lengths  so  that  Eqn.(80)  could  be  used  to 
determine  e (v) . That  is,  a regression  analysis  was  used  to 
obtain  the  slope  of  the  in  T vs  b data.  The  results  are 
presented  in  Figs.  5 a and  b where  it  can  be  seen  that  the 
value  of  e(v)  is  a function  of  concentration  in  low  absorp- 
tion regions,  but  is  essentially  independent  of  concentration 
for  e(v)  > - 10*  liters/mole-cm,  i.e.  at  high  absorption. 

It  should  be  remembered  here  that  spectrophotometric  error  is 
greatest  in  regions  where  transmission  is  either  very  small 
(i.e.  T < - 10%)  or  very  large  (T  > ~ 90%).  The  data  re- 
flecting the  concentration-dependence  corresponds  to  high 
transmission  values  whereas  the  peak  concentration-indepen- 
dent values  of  e(v)  were  obtained  from  nearly  optimum  trans- 
missions. Thus  for  e(v)  > ~ 10*,  the  results  have  a higher 
confidence  factor  than  for  lesser  absorption  regions. 
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Preliminary  measurements  at  TBP  concentrations  higher 
than  those  reported  in  Fig.  5 a and  b have  been  made.  They 
indicate  that,  for  the  very  strong  band  near  2950  cm  the 
value  of  e(v)  sharply  drops  at  c ~ 0.3M  by  a factor  of  ap- 
proximately 1/3  and  remains  essentially  constant  out  to  the 
concentration  of  the  pure  liquid  (i.e.  c = 3.65M).  Measure- 
ments aimed  at  corroborating  this  effect  as  well  as  pro- 
viding quantitative  data  on  the  pure  liquid,  will  be  reported 
in  the  future. 


C.  Cellulose 

Samples  of  regenerated  cellulose  film  were  obtained  from 
the  Film  and  Packaging  Division  of  FMC  Corporation.  Two  film 
thicknesses  were  specially  prepared  for  our  use  and  contained 
no  plasticizer.  The  films  were  dried  under  tension  in  FMC 
Corporation's  Laboratory.  While  the  moisture  content  of  cel- 
lulose is  very  sensitive  to  atmosphere  humidity,  we  were  un- 
able to  determine  any  appreciable  difference  between  the 
spectrum  of  a film  stored  in  air  and  one  which  had  been  dryed 
for  24  hours  under  vacuum.  No  additional  drying  at  elevated 
temperatures  was  attempted  by  us. 

Cellulose  is  a naturally  occurring  material  having  a typi- 
cal fibre  structure.  It  is  the  main  component  of  the  cell 
walls  of  most  plants  and,  thus,  has  been  estimated  to  comprise 
about  33%  of  all  vegetable  matter.  Cellulose,  together  with 
starch  and  glycogen  are  polysaccharides  composed  of  glucose 
only.  The  glucose  units  are  linked  by  6-glycoside  bonds  re- 
sulting in  a repeat  or  monomer  unit  represented  by 
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and  referred  to  as  a cellobiose  unit. 

Industrial  cellulose,  as  used  here,  is  extracted  from 
wood  pulp  by  a straight-forward,  but  carefully  controlled  pro- 
cess. Our  laboratory  is  not  equipped  to  perform  the  extrac- 
tion procedure.  The  films  kindly  prepared  for  us  by  FMC 
Corporation  were  received  with  average  thicknesses  of  0.45 
mil.  (.0011  cm)  and  0.95  mil.  (.0024  cm).  Normally,  such 
optical  path  lengths  for  transmission  measurements  are  suit- 
able for  determining  absorption  coefficients.  However,  cel- 
lulose is  highly  absorbing  in  the  3600-3000  cm”*  (2.8-3.33vj) 
and  1200-950  cm"*  (8.3-10.5p)  regions.  This  is,  of  course, 
a highly  desirable  property  in  terms  of  the  goals  of  the 
overall  project.  The  low  transmission  of  the  samples  in  these 
regions  made  quantitative  measurements  impractical.  Therefore, 
in  Figs.  6 a and  6 b,  where  the  behavior  of  a(v)  for 
cellulose  is  presented,  those  regions  for  which  numerical 
values  are  not  reported  are  indicated  by  dashed  lines.  It 
may  be  possible  to  produce  samples  of  cellulose  that  are  thin 
enough  to  yield  transmissions  of  30%  in  those  highly  absorb- 
ing regions,  yet  are  mechanically  stable  enough  to  be  treated 
as  free  films.  However,  to  do  so  should  require  a special 
order  to  FMC  Corporation  without  additionally  imposing  on 
their  generosity  and  cooperation. 

D.  Polyvinyl  Alcohol 

Polyvinyl  alcohol  (PVA)  resin,  in  powdered  form,  was  do- 
nated by  American  Hoechst  Corp.  Material  received  was  de- 
signated as  Moviol  4-88  and  characterized  by  a viscosity  of 
4 cP  (centipoise) , for  a 4%  aqueous  solution  at  20®C, and  an 
88%  degree  of  hydrolysis.  PVA  solutions  were  produced  by 
adding  the  powder  to  stirred  cold  water,  then  heating  the 
mixture  to  85-90®C  until  dissolution  was  complete.  The  re- 
sulting solution  is  then  poured  into  casting  dishes  upon 
which  films  are  formed  when  the  solvent  has  evaporated. 
Approximately  1%  solutions  (by  weight)  were  used.  Sample 
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thickness  was  controlled  by  varying  the  volume  of  solution 
poured  into  the  casting  dish. 

Three  films  were  produced  with  thicknesses  of  23.6±2.5ym, 
3.81i0.00um,  and  2.90±0.61ijm.  The  errors  in  the  thicknesses 
were  determined  by  measuring  only  the  variations  of  those 
areas  of  the  films  "seen"  by  the  optical  beam  of  the  in- 
strument. The  optical  density  of  the  23. Gym  film  was  too 
large  to  allow  for  quantitative  results  to  be  reported  over 
the  wavelength  ranges  of  interest.  Therefore,  values  of 
a(v)  reported  in  Fig. 7 a and  b represent  data  obtained  from 
ttie  two  thinner  films. 

E.  Ureaformaldehyde 

Powdered  ureaformaldehyde  (UF)  resin  was  supplied  by 
R,  W.  Doherty  (CSL-Edgewood  Arsenal)  and  characterized  as 
having  an  approximate  particle  size  of  4ym  diameter.  The 
resin  is  crosslinked*  and  so  is  not  soluble.  Because  of  its 
relatively  poor  characterization,  ureaformaldehyde  could  only 
be  investigated  as  a powder  suspended  in  a matrix  (Nujol  Oil) . 
Quantitative  measurement  techniques  for  such  powders  are 
presently  being  developed  in  our  laboratory.  Consequently, 
only  a qualitative  spectrum  of  UF  is  presented  (see  Fig. (8a  and  8b) 
at  this  time. 


♦Personal  communication  with  J.  Schwendeman,  Monsanto  Corp. 
Dayton,  Ohio 
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Fig.  6b  - Absorption  Coefficient 
of  Cellulose.  6-13um. 


51 


a (v) , cm 


-Cl 


Ureaformaldehyde . Sample 
prepared  as  a KBr  pellet. 


GLOSSARY  OF  SELECTED  TERMS 


Symbol 

Quantity 

Units 

Definition 

A 

Wavelength 

Viir  = micrometers 

V 

Wavenumber 

-1 

cm 

c 

II 

>'|i-‘ 

c 

Speed  of  light 

meters 

sec 

c = i//T~u  = 

0 0 

3x10®  meters/s 

f 

Frequency 

Hertz,  Hz 

t-tl 

II 

(1) 

Angular  Frequency 

radians 

sec 

(1)  = 2Trf  = 2itCv 

K 

Wave  Vector, 
Magnitude 

radians 

cm 

_ 271 

^ ~ X 

K*  (uj) 

Wave  Vector, 
(Complex) 

radians 

cm 

K*  = ^ n* 

n 

Index  of 

Refraction 

none 

_cfv  = speed  of 
v[  light  in  medium 

n*  (V) 

Index  of  Refrac- 
tion (Complex) 

none 

n* (v)  = n (v) +ik (v) 

k(v) 

Extinction 

Coefficient 

none 

k (v)  = Im[n* (v) ] 

A 

Absorbance  or 
Optical  Density 

none 

^ ” ^°^10  ~f~ 

T 

Transmission 

none 

I /ratio  of  \ 

I 1 transmitted! 
o|  to  incident! 
y power  / 

a (v) 

Absorption 

Coefficient 

cm  ^ 

a (v)=gfa  (v)=^^  (v) 

^ \ 

= 477vk  (v)  1 

Z 


Symbol 

Quantity 

Units 

b 

Optical  Path 
Length 

cm 

e(v) 

Molar 

Extinction 

Coefficient 

liter 

mole-cm 

Concentration 

mole 

c 

liter 

Definition 


£{v) 


A 

be 


|e  (v)  = (v) 
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Appendix  A:  Development  of  General  Wave  Equation  to  Describe 
The  Propagation  of  Light  in  Dielectric  or 
Conductive  Media. 


The  behavior  of  electromagnetic  waves  in  any  medium  is  des- 
cribed by  the  classical  wave  equation,  which  in  turn  is  derived 
from  Maxwell's  Equations.  In  this  section  the  development  of 
the  appropriate  wave  equation,  of  interest  to  this  work,  will 
be  briefly  presented. 

The  material  parameters  used  to  characterize  the  medium 
with  respect  to  electromagnetic  phenomena  are,  in  SI  units,  the 
relative  permitivity  = e/€^le^=  8.85  x 10”'^  farad  m~‘  = 

permitivity  of  free  space)  and  the  relative  permeability 
= ij/vj^(u^=  1.26  X 10”®  henry  m”*  = permeability  of  free 
space).  The  velocity  of  propagation  of  the  wave  in  this  medium 
is 


V 


/k  e /c  u 
e 0 mo 


(A-1) 


which,  for  free  space  (in  vacuo  K = x =1)  reduces  to 

e m 


= 3 X 10®  m/s 


The  index  of  refraction  of  the  medium  is  given  by 


(A-2) 


u = -g-  = (A-3) 

V e m 

For  most  transparent  optical  media  k = 1 so  that 

m 

This  result  is  generally  only  valid  for  dilute  gases  and  non- 
polar solids,  where  high  static  polarizability  is  present.  In 
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addition,  n is  a function  of  frequency  (dispersion)  so  that  to 
characterize  the  medium  optically,  these  effects  must  be  account- 
ed for  and  (A-4)  is  seen  to  be  a special  limiting  case. 

Maxwell's  Equations  for  a non-magnetic , electrically  neu- 
tral medium  can  be  reduced  to 


V X 2 = - 


7 X H = e + -1^  + D 

0 3 t 3 t 


(A-5a) 


(A- 5b} 


V . E = - 


(A-5c) 


7 • H = 0 


(A-5d) 


where  E is  the  electric  field  intensity,  H is  the  magnetic  field 
intensity, 

? = (e-e^)  E = (A-6) 


is  the  electric  polarization  (y 
bility)  and 

-t  -t- 

3 = OE 


- 1 = electric  suscepta- 


(A-7) 


is  the  current  density  (o  ® electrical  conductivity) . 

By  manipulating  the  ad>ove  equations  [take  the  curl  of  (A-5a) 
and  time  derivative  of  (A-5b)  and  eliminate  H between  the  two 
resulting  equations]  one  can  get 

7 X (7xi)  + =■  - y - y^  ^ (A-8) 


C*  3 t* 


° 3 t* 


From  the  vector  identity 


7 X (7  X E)  - 7(7*l)  - 7*  2 


(A-9) 


and  the  fact  that  V*E  = 0 for  an  electrically  neutral  media,  it 
follows  that 


V^E 


1 

1 t^ 


+ 


ll 

3 t 


(A-10) 


The  two  terms  on  the  right  hand  side  of  (A-10)  are  called 
source  terms  and  are  interpreted  as  follows: 


(a)  For  a non-conducting  medium,  accounts  for 

3 t ^ 

the  effect  of  polarization  charges  and  is  the 

important  source  term,  being  responsible  for  optical 
effects  such  as  dispersion,  absorption,  etc. 

3 1 

(b)  For  a conducting  medium,  is  the  important  term 

and  is  responsible  for  the  high  reflectivity  and  optical 
opacity  of  metals. 

(c)  For  semiconductors,  both  source  terms  must  be  considered. 

Equation  (A-10)  is  the  starting  point  for  discussions  of  dis- 
persion and  metallic  absorption  considered  in  this  work. 


